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Using first-principles calculations within density functional theory, we investigate the intrinsic 
spin Hall effect in monolayers of group- VI transition-metal dichalcogenides MX2 (M = Mo, W and 
X = S, Se). MX2 monolayers are direct band-gap semiconductors with two degenerate valleys 
located at the corners of the hexagonal Brillouin zone. Because of the inversion symmetry breaking 
and the strong spin-orbit coupling, charge carriers in opposite valleys carry opposite Berry curvature 
and spin moment, giving rise to both a valley-Hall and a spin-Hall effect. We also show that the 
intrinsic spin Hall conductivity in inversion-symmetric bulk dichalcogenides is an order of magnitude 
smaller compared to monolayers. Our result demonstrates monolayer dichalcogenides as an ideal 
platform for the integration of valleytronics and spintronics. 



PACS numbers: 72.25.Dc,73.63.-b,75.70.Tj,78.67.-n 

I. INTRODUCTION 

In crystalline semiconductors, it often happens that 
the conduction band minima and valence band maxima 
are located at degenerate but inequivalent valleys. Well- 
known examples include graphene, 1 bismuth thin films, 2 
and AlAs quantum wells. 3 Since the valleys are usually 
separated by a large distance in momentum space, inter- 
valley scattering is greatly suppressed in the presence of 
smooth scattering potential, rendering the valley index 
an intrinsic property for low-energy carriers. Motivated 
by this observation, there has been a growing interest 
in exploiting the valley index in electronic devices, much 
in the same way as the spin index is used in spintronic 
applications. This is the subject of valleytronics. 

Recently, a general scheme based on inversion symme- 
try breaking has been proposed to generate and manip- 
ulate the valley polarization. 4,5 The central idea is that 
under inversion symmetry breaking, the valley index can 
be associated with distinctive physical quantities such as 
the Berry curvature and orbital magnetic moment. 6 Us- 
ing graphene as an example, the authors of Ref. 4 and 5 
showed that inversion symmetry breaking allows a valley 
Hall effect in which carriers in different valleys flow to op- 
posite transverse edges when an electric field is applied, 
leading to a finite valley polarization along the edges. 4 
Furthermore, it also gives rise to valley-contrasting cir- 
cular dichroism in the momentum space, which takes the 
extreme form of optical selection rules at high symmetry 
points. 5 Other approaches have also been proposed. 1 ' 7 
However, they all rely on carefully prepared geometry at 
the atomic scale, which is difficult to control in experi- 
ments. 



In general, inversion symmetry breaking also lifts the 
spin degeneracy of energy bands in the presence of spin- 
orbit coupling (SOC). As required by time-reversal sym- 
metry, the spin-splitting in opposite valleys must be op- 
posite, therefore the valley carriers can be also distin- 
guished by their spin moments. This is the basis of cou- 
pled spin and valley physics. However, the SOC is negligi- 
bly small in graphene, 8,9 preventing further investigation 
along this direction. In a recent work, 10 we have studied 
monolayers of group- VI transition-metal dichalcogenides 
for the following reasons: (i) the inversion symmetry is 
explicitly broken in monolayers; (ii) the conduction and 
valence bands of these materials harbors a multi-valleyed 
structure; 11,12 and (iii) the SOC is substantial due to the 
presence of heavy metal atoms. 13 Therefore these materi- 
als provide a perfect platform to investigate the interplay 
between spin and valley degrees of freedom. Based on an 
effective k • p model, we predicted that the valley Hall 
effect is accompanied by a spin Hall effect in both n- 
and p-doped systems, and the valley-dependent optical 
selection rule also becomes spin-dependent. 10 

Monolayers of group- VI dichalcogenides also display 
excellent optical properties for practical applications. 
Recent experiments have demonstrated that M0S2, a pro- 
totypical group- VI dichalcogenide, crossovers from an in- 
directgap semiconductor at multilayers to a direct band- 
gap one at monolayer. 14,15 The direct band-gap is in 
the visible frequency range, most favorable for optoelec- 
tronic applications. Experimental evidence of the valley- 
dependent optical selection rule in monolayer M0S2 has 
been recently reported based on polarization-sensitive 
photoluminescence measurement . 16 18 

In this work we present a comprehensive first-principles 
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study of the coupled spin and valley physics, focusing on 
the Hall effects of valley and spin. We show that, be- 
cause of the inversion symmetry breaking and the strong 
SOC, charge carriers in opposite valleys carry opposite 
Berry curvature and spin moment, giving rise to both the 
valley- and spin-Hall effect. Our first-principles calcula- 
tions provide a quantitative basis for the k ■ p model de- 
rived in Ref. 10, and subtle differences between these two 
are discussed. We also show that the intrinsic spin Hall 
conductivity in inversion-symmetric bulk dichalcogenides 
is an order of magnitude smaller compared to monolay- 
ers. Our result demonstrates monolayer dichalcogenides 
as an ideal platform for the integration of valleytronics 
and spintronics. 



II. METHODOLOGY 

The electronic ground-state calculations in this work 
were performed using full-potential linearized aug- 
mented plane- wave method, 19 implemented in the pack- 
age wien2k. 20 Exchange-correlation effect was treated 
with the Perdew, Burke, and Ernzerhof parameterized 
generalized-gradient approximation. 21 The crystal struc- 
ture were adopted from the first-principles optimized re- 
sults. 13 For the slab model, a 20 A thick vacuum layer was 
used to avoid the interactions between adjacent mono- 
layers. The converged ground-states were obtained using 
/c-mesh 16 x 16 x 3 for bulk and 16 x 16 x 1 for monolayer 
in the first Brillouin zone, both with i^max^MT = 7.0, 
where Rmt represents the smallest muffin-tin radius and 
K max is the maximum size of reciprocal-lattice vectors. 
Wave functions and potentials inside the atomic sphere 
were expanded in spherical harmonics up to /=10 and 
4, respectively. Spin-orbit coupling was included by a 
second- variational procedure, 19 where states up to 9 Ry 
above Fermi level were included in the basis expansion. 

To calculate the Berry curvature and Hall conductiv- 
ity, we first computed the Wannier functions by the max- 
imally localized algorithm, 22,23 implemented in the pack- 
age WANNIER90. 24 The construction of maximally local- 
ized Wannier functions is a non-self-consistent process 
on a uniform 8x8x8 grid of /c-point with formerly 
converged self-consistent charge potential. Transition- 
metal dichalcogenides have the chemical formula MX^ 
(M = Mo, W and X = S, Se). In the case of monolay- 
ers, there are 22 bands in the energy range from about -6 
to 5 eV, mainly formed by M d- and X p-orbitals. With 
this in mind, we chose ten d orbitals on atom M and six p 
orbitals on each atom X as the initial guess of the Wan- 
nier functions. After less than 200 iterative steps, the 
total Wannier spread was well converged down to 10 _T 
Bohr 2 . On the other hand, in bulk systems, the unit 
cell contains two formula units, and there are 44 bands 
also formed by M d- and X p-orbitals. The construction 
process in bulk is similar to that in monolayers except 
for WS2 and WSe2- In these materials, the conduction 
bands are entangled with higher bands. The disentan- 
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FIG. 1. (Color online) (a) Side view of the unit cell of the 2H- 
MX 2 structure with M=Mo, W and X=S, Se. It contains 
two MX2 monolayers separated by a Van der Waals gap. (b) 
Top view of the MX2 monolayer. The black lines indicate 
the unit cell in ab plane, (c) The first Brillouin zone and high 
symmetry points of the MX2 monolayer, (d) The Wannier 
functions of the MX2 monolayer, including five d-orbitals on 
M atom and three p-orbitals on each X atom. 



glement approach 23 was applied to bulk WS2 and WSe2- 
Once the Wannier functions were obtained, we followed 
Ref. 25 to calculate the Berry curvature and integrate it 
over the Brillouin zone to obtain the Hall conductivity. 



III. BAND STRUCTURE AND BERRY 
CURVATURE 

In this section we present the electronic band structure 
and Berry curvature of monolayers of group- VI dichalco- 
genides, using M0S2 as an example. We show that inver- 
sion symmetry breaking gives rise to two physical quanti- 
ties, the spin moment and the Berry curvature, that can 
be used to distinguish valley carriers. 

Structurally, M0S2 can be regarded as strongly bonded 
two-dimensional S-Mo-S layers that are loosely coupled 
to one another by Van der Waals interactions. Within 
each layer, the Mo and S atoms form hexagonal lattices 
in separate planes with each Mo atom coordinated by six 
nearest-neighboring S atoms in the trigonal prismatic ge- 
ometry [Fig. 1(a) and (b)]. In its bulk form, M0S2 has the 
so-called 2H stacking order with space group P6s/mmc 
(Dg^), which is inversion symmetric [Fig. 1(a)]. Be- 
cause of the weak interlayer coupling, this layered com- 
pound can be easily exfoliated into monolayers by me- 
chanical 26,27 and chemical 28 means. In monolayer M0S2, 
the space group is reduced to P6ra2 (D\ h ) with explicit 
breaking of inversion symmetry. 

Before moving on to the discussion of coupled spin and 
valley physics, we briefly study the orbital characters via 
the Wannier functions. The partial density of states of 
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FIG. 2. (Color online) The calculated electronic structure of 
M0S2 with the spin-orbit coupling, (a) The partial density of 
states for Mo-4d and S-3p orbit als, respectively, in the unit of 
states/eV/cell. (b) The band structure with the projection of 
spin operator s z (color map). The red and blue colors indi- 
cate the spin-up and -down states, respectively. The optical 
transitions between the VBM and the CBM are coupled exclu- 
sively with cr+ (a—) circular polarizations at the inequivalent 
valleys K (K f ). 10 



monolayer M0S2 shown in Fig. 2(a) clearly shows that 
in the energy range of -6 to 5 eV the contributions to 
electronic states mainly come from the Mo d- and S p- 
orbitals, whereas other orbitals have vanishing contribu- 
tions. In Fig. 1(d), we plot the typical Wannier functions 
for five Mo d- and three S p-orbitals together with the 
unit cell in real space. If the SOC is turned on, the num- 
ber of orbitals will be doubled and altogether 22 bands 
are formed. Under the crystal field of trigonal prismatic 
coordination, the d-orbitals split into {^2}, {d xz ,d yz }, 
{d x 2_ y 2,d xy } and p-orbitals split into {p z }, {Px,Py}- At 
the two inequivalent valleys, K and K' \ the valence band 
maximum (VBM) is constructed by the Mo {d x i_ y i,d xy ) 
orbitals with some mixing from the S {p x ->Py} orbitals, 
while the conduction band minimum (CBM) is domi- 
nated by Mo d z 2 orbitals. These orbital characters of 
band-edges are in line with the analysis of the k -p model 
in Ref. 10. 

Figure 2(b) shows the fully relativistic band structure 
of monolayer M0S2 with the projection of spin operator 
s z , i.e., (ip n k \ s z |^nfc) 7 obtained from Wannier interpola- 
tion. 29 We can see that there is a direct band-gap at the 
two inequivalent corners K and K' of the Brillouin zone. 
Furthermore, a large spin splitting (~ 150 meV) appears 
at the VBM with opposite spin moments at the two val- 
leys, as a result of inversion symmetry breaking. 13 This 
indicates that in addition to their valley index, the valley 
carriers in the valence bands can be also distinguished by 
their spin index. On the other hand, since the CBM state 
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FIG. 3. (Color online) The Berry curvatures of monolayer 
M0S2 along the high symmetry lines (a) and in the 2D /c-plane 
(b). The spin Berry curvatures of monolayer M0S2 along the 
high symmetry lines (c) and in the 2D fc-plane (d). All of the 
Berry curvatures are in the atomic unit (Bohr 2 ). 



is made of the Mo d z i orbital, SOC is inactive and the 
CBM remains degenerate (to first order of the SOC). 

In the presence of inversion symmetry breaking, the 
charge carriers also acquire a valley-contrasting Berry 
curvature. 4,6 According to the Kubo-formula, 30 ' 31 the 
Berry curvature Q (k) of the occupied states can be writ- 
ten as 



n(k) = Y,f^ (fe), 



(i) 



and 



.(*) = - E 



21m (l/>nk\Vx iVVfc) {^n'k\Vy \^nk) 

(E n > — E n ) 2 



(2) 

where \ip n k) is the Bloch function with the eigenvalue E n , 
f n the Fermi-Dirac distribution function, and v x ^ y \ the 
velocity operators. Here, we have used the maximally 
localized Wannier functions as the basis to calculate the 
Berry curvature [Eq. (2)] and the spin Berry curvature 
[below in Eq. (6)]. 25 Figure 3(a) shows fl(k) of mono- 
layer M0S2 along the high-symmetry lines. We can see 
that fi(fc) is significantly peaked at both K and K' but 
with opposite signs. The /c-space contrasting fi(fc) in 
systems without inversion symmetry is a key quantity to 
characterize the chirality of the Bloch electrons and is the 
basis for valley-contrasting phenomena. 4,5,10 Away from 
the two valleys, fi(fe) decays rapidly and vanishes at the T 
and M points. We also plot the map distribution of Q(k) 
in the 2D /c-plane, as shown in Fig. 3(b), which clearly 
shows the C3 symmetry of the system. This is in con- 
trast with an energy counter plot, which would display 
the Cq symmetry, i.e., the two valleys are energetically 
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curvatures. Hence, the total anomalous Hall conductiv- 
ity vanishes because of time reversal symmetry. If a fi- 
nite valley polarization can be generated, for example, by 
shining the sample with circularly polarized light, then 
a charge Hall current will appear. 10 On the other hand, 
since both the valley and spin current remains invariant 
under time-reversal, the valley Hall and the spin Hall 
effect can appear in time-reversal invariant systems, as 
long as the inversion symmetry is broken. 32 

Finally, we compare the value of the Berry curvature at 
the VBM and CBM from both the effective k • p model 10 
and first-principles calculation in Table I. The excellent 
agreement between them further confirms the validity of 
the k • p model. 



IV. THE INTRINSIC SPIN HALL EFFECT 



FIG. 4. (Color online) (a) The intrinsic spin Hall conductivity 
^xy ( e V^) as a function of the Fermi energy for monolayer 
M0S2. The energy zero point (true Fermi level) is at the 
VBM. Two large peaks close to the CBM and VBM are de- 
noted by Pi and P 2 , respectively. (b)(c) The band structure 
(up panel) and the spin Berry curvature (down panel) when 
the Fermi level shifts to the positions of the peak Pi and P 2 , 
respectively. There is a small peak of Q s with negative value 
along the M-Y line in (c). (d) The n- and p-doping charge 
as a function of the Fermi energy. The low doping regimes 
just above the CBM and below the VBM (indicated by red 
circles) are more relevant in experiments. 



indistinguishable, but they can be distinguished by their 
Berry curvatures. 

The Berry curvature drives an anomalous transverse 
velocity in the presence of an electric field E: 6 



h 



e x n(fe) , 



(3) 



which is responsible for the intrinsic contribution to the 
anomalous Hall effect. 31 However, in our systems, the 
charge carriers in the two valleys have opposite trans- 
verse velocities due to the opposite signs of the Berry 



TABLE I. Comparison of the Berry curvatures at VBM and 
CBM between the present first-principles calculations (the 
first line) and the k • p model in Ref. 10 (the second line). 
Qv(c)t(l) i s the Berry curvature of the valence (conduction) 
band with spin t (4<) 5 gi ven i n the unit of Bohr 2 . 





M0S2 


MoSe 2 


WS 2 


WSe 2 




38.8 


39.7 


59.8 


64.3 
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36.5 


55.4 


60.0 




31.6 


30.0 


34.9 


34.7 




29.5 


28.4 


34.2 


33.3 




-35.7 


-36.8 


-54.7 


-59.2 




-35.3 


-36.5 


-55.4 


-60.0 




-28.8 


-27.3 


-31.0 


-30.8 




-29.5 


-28.4 


-34.2 


-33.3 



As discussed above, both the valley Hall and spin 
Hall effect exist in MX2 monolayers due to the valley- 
contrasting Berry curvature. Note that the valley index 
is defined only in the vicinity of the valleys, whereas the 
spin index is defined everywhere in the Brillouin zone. 
Therefore we will only calculate the intrinsic spin Hall 
conductivity. For hole-doped samples, when the Fermi 
energy lies between the spin-split VBM states, the valley 
Hall conductivity coincides with the spin Hall conductiv- 
i.v. 1 " 

At zero-temperature and clean limit, the intrinsic spin 
Hall conductivity (ISHC) tensor is given by 



Va (2tt)' 



(4) 



For the convenience of discussion, in the following, we 
multiply a factor 2e/hto the calculated ISHC to convert 
its unit to charge conductivity. We can carry out the 
calculation of a x again using the Kubo-formulas 33 ' 34 



n s (fc) = 2/„n-(fc), 



(5) 



and 



n s (k) = 2Im i^nkljx h/Vfc) {^n'k\Vy \^nk) ^ 



(E n t — E n ) 



where j x is the spin current operator defined as \ {s z v x + 
v x$z)- We add a superscript s for the spin Berry curva- 
ture in order to distinguish them from the ordinary Berry 
curvature in Eq. (1) and (2). We can see that Q s (k) of 
monolayer M0S2 is peaked at both K and K' with the 
same sign, as shown in Fig. 3(c). This can be under- 
stood as the following. At the VBM, s z remains a good 
quantum number, and the spin Berry curvature is simply 
given by n s (K) = s z Q(K). As both s z and Q(K) flip 
sign when K — )► K . Q S (K) remains the same. Fig- 
ure 3(d) shows the map distribution of Q s (k) in the 2D 
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FIG. 5. (Color online) The intrinsic spin Hall conductivity 
a xy in the low doping regimes for the monolayer (a)-(d) and 
bulk (e)-(h) of the MoS 2 , MoSe 2 , WS 2 , and WSe 2 , respec- 
tively. Dashed lines indicate the band-edges. Note that the 
band-gaps of bulk are smaller than the monolayer ones. The 
unit of a xy is e 2 /h (= 2.43 x 10 _4 n _1 ) for 2D monolayer sys- 
tem, whereas Q -1 cm -1 for 3D bulk system. For quantita- 
tively comparing the a xy in bulk and monolayer, one needs to 
divide the a xy in monolayer by its thickness. 



fc-plane. We observe that it has a clear Cq symmetry 
rather than the C3 symmetry in Fig. 3(b). 

By integrating Q s (k) over the occupied states, we ob- 
tain G xy as a function of Fermi level for monolayer M0S2, 
shown in Fig. 4(a). Here, we set the energy zero point 
(true Fermi level) at the VBM, and calculate G xy by 
rigidly shifting the Fermi level position. For n-doped 
monolayer M0S2, the calculated o s sharply reaches its 
maximum value of 0.89 e 2 /h at 2.19 eV. When further 
increasing the doping concentration, G s xy displays a com- 
plex behavior with both dramatic oscillations and sign 
changes, but it eventually goes to zero above 4.5 eV. For 
p-doped monolayer M0S2, o 8 xy has two large peaks, re- 
spectively, with positive value of 0.29 e 2 /h at —1.52 eV 
and negative value of —0.32 e 2 /h at —2.64 eV. In order to 



TABLE II. The ISHC a s xy of monolayer MX 2 calculated at 
hole-doping concentration n/ l =1.0xl0 13 cm -2 . The slopes 
of a xy when Fermi level lies inside the spin splitting gaps 
below the VBM are also listed. The first and second lines are 
obtained from first-principles calculation and two-bands k • p 
model 10 , respectively. 





MoS 2 


MoSe 2 


ws 2 


WSe 2 


slope (e 2 /h/eV) 


0.21 


0.24 


0.17 


0.18 




0.20 


0.23 


0.20 


0.23 


a s xy (10- 2 e 2 /K) 


-0.58 


-0.92 


-1.13 


-1.30 


-0.57 


-1.00 


-1.33 


-1.52 



analyze the cause of the large peaks, we take Pi and P2 
as examples, indicated in Fig. 4(a). When the Fermi level 
is located at the positions of P\ and P2, the band struc- 
tures and spin Berry curvatures are shown in Fig. 4(b) 
and 4(c), respectively. We can see that Q s (k) is often 
peaked at the places where the Fermi level crosses some 
tiny band-gaps induced by the SOC. Similar behavior of 
Q s (k) for the other large peaks of the G xy are also found. 
The appearance of large peaks of the positive (negative) 
Q s (k) leads to the positive (negative) peaks of the g s . 

Although giant G xy can be realized at those peak po- 
sitions, such a high level of doping is unrealistic in ex- 
perimental situations. As indicated in Fig. 4(d), the 
Fermi level position of Pi and P2 are 2.19 eV and —1.52 
eV, respectively, corresponding to electron concentration 
n e =1.73 e/cell (= 1.96xl0 15 cm -2 ) and hole concentra- 
tion n /l =3.18 e/cell (^ 3.60xl0 15 cm" 2 ). This could 
be difficult in experimental conditions either by chemical 
adsorption or by gate voltage. For example, the high- 
est carrier concentration in two-dimensional graphene is 
only up to 10 13 cm -2 (Ref. 26 and 35). In contrast, the 
low doping regimes just above the CBM and below the 
VBM are more relevant in experiments, as indicated by 
red circles in Fig 4(a) and 4(d). In the following, we 
only focus on this regime. 

Figure 5(a)-(d) show the G s xy in the low doping regimes 
for monolayer M0S2, MoSe2, WS2, and WSe2, respec- 
tively. We can clearly see that G xy increases with the 
SOC strength as the atoms becomes heavier, and the p- 
doped samples generally has larger G s xy than the n-doped 
samples. This is due to the large spin splitting at the 
VBM [see Fig. 2(b)]. Focusing on the valence bands, we 
find G xy is proportional to the Fermi energy. The ex- 
tracted slopes are listed in Table II together with the 
results from the k -p model. 10 . When the Fermi level lies 
inside the spin splitting gaps, the hole concentration is on 
the order of 10 13 cm -2 , which is realistic for experiments. 
Taken rih=l. Ox 10 13 cm -2 as an example, we list the cal- 
culated G xy from both present first-principles calculation 
and the k-p model 10 in Table II. The ISHC G s xy listed 
in Table II are also comparable to those in p(n)-doped 
semiconductors GaAs, Si, Ge, and AlAs. 33,34 

Two remarks are in order. First, here we only com- 
pared the ISHC for p-doped samples for both first- 
principles calculations and the k-p method. For n-doped 
samples, the situation is more complicated. As we can 
see in Fig. 6, in both WS2 and WSe2 the conduction band 
has a second local minimum between T and K, which is 
very close to the band-edge at K. Even under light dop- 
ing (~ 1.0 x 10 13 cm -2 ), both minima will be occupied 
and contribute to the total ISHC, rendering the compar- 
ison between first-principles and k • p method meaning- 
less. This shows the limitation of the k • p method: it 
only captures the physics around the K point, and first- 
principles study give us a more complete picture. We also 
note that there is a small spin splitting at the CBM for 
both WS2 and WSe2- This is due to the much heavier 
W atom compared to Mo (SOC scales as Z 4 , where Z is 
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FIG. 6. (Color online) The band structure of WS 2 and WSe 2 
monolayers. The CBM of WSe 2 monolayer still locates at K 
point. The dashed line indicate the position of the Fermi level 
at carrier concentration of 1.0 x 10 13 cm -2 for both p- and 
n-doped samples. 



the atomic number). In the k-p model, this splitting can 
be taken into account by considering the second-order ef- 
fect due to the SOC. Second, one may notice that a xy is 
nonzero in the band-gap. These nonzero values are not 
due to numerical errors, but actually reflect the finite hy- 
bridizations in real materials, similar to what has been 
reported in GaAs and Si. 33 The MX2 monolayer studied 
here can be viewed as a generalization of the concept of 
spin Hall insulator proposed by Murakami et a/., 36 such 
as PbTe, which is a conventional band insulator but has 
nonzero a xy without any doping. 

As a comparison we also calculate a xy for bulk M0S2, 
MoSe 2 , WS 2 , and WSe 2 , shown in Fig. 5(e)-(h). The 
shapes of cr xy for all bulk dichalcogenides look rather sim- 
ilar, though the magnitudes and positions of the peaks 
may differ. In bulk, a xy is zero immediately below the 
VBM, which is different from the immediate increasing 
of the <J xy in monolayers. The reason is that the VBM in 
bulk always locates at the Y point, which has negligibly 
small contribution to a xy . To compare with monolayers, 
we divide the o s in monolayer by its thickness and find 
that at the same Fermi level the a xy in bulk is about an 
order of magnitude smaller than that in monolayers. 

Finally we mention that here we only considered the 



intrinsic contribution to the spin Hall effect, in which the 
spin Hall current is driven by the Berry curvature of the 
Bloch bands. There are also extrinsic contributions com- 
ing from scattering of impurities and phonons. When the 
sample is hole-doped, the effect of phonon scattering on 
the SHC should be weak because the phonon scattering 
will mostly contribute to intra- valley scattering, in which 
the spin z-component is nearly conserved due to the large 
spin-orbit splitting at the valence band top. On the other 
hand, impurity scattering can provide the large momen- 
tum transfer needed for the inter-valley scattering, and 
their effect on valley-dependent transport properties re- 
mains to be investigated. 

V. SUMMARY 

In summary, using the first-principles calculations, we 
have investigated the intrinsic spin Hall effect in mono- 
layers M0S2, MoSe2, WS2, as well as WSe2, driven by 
valley-contrasting Berry curvature. We find that the 
ISHC is comparable to that in p(n)-doped semiconduc- 
tors GaAs, Si, Ge, and AlAs. 33,34 We show that the ef- 
fective model may not be adequate to describe the low- 
energy dynamics in WS2 and WSe2 monolayers. We also 
calculated the ISHC in inversion-symmetric bulk systems 
and find that it is about an order of magnitude smaller 
than the ISHC in monolayers. The large ISHC, plus other 
interesting physical properties of these materials, such as 
giant spin splitting at VBM 13 and valley- selective circu- 
lar dichroism, 10,16 18 characterize these materials as an 
exciting platform for the application of the valleytronics 
and spintronics. 
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